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Abstract 


Differential cross sections as a function of transverse momentum pj are presented 
for the production of Y(nS) (n = 1, 2, 3) states decaying into a pair of muons. Data 
corresponding to an integrated luminosity of 4.9 fb~^ in pp collisions at = 7 TeV 
were collected with the CMS detector at the LHC. The analysis selects events with 
dimuon rapidity \y\ < 1.2 and dimuon transverse momentum in the range 10 < 
Pi < 100 GeV. The measurements show a transition from an exponential to a power- 
law behavior at pj ~ 20 GeV for the three Y states. Above that transition, the Y(3S) 
spectrum is significantly harder than that of the Y(1S). The ratios of the Y(3S) and 
Y(2S) differential cross sections to the Y(1S) cross section show a rise as pj increases 
at low Pj, then become flatter at higher pj. 
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1 Introduction 

Hadronic production of S-wave bb mesons has been extensively studied for many years. At the 
CERN LHC, the CMS HHIH, ATLAS HSj, and LHCb |4| Collaborations have published results on 
Y (nS) (n = 1, 2, 3) production cross sections times dimuon branching fractions in pp collisions 
at \/s = 7TeV as a function of the Y transverse momentum pj, rapidity y, and polarization [Sf. 
The CMS and ATLAS px and \y \ distributions in the central rapidity region \y\ < 2.0 are similar 
in shape to those from pp production at ^/s = 1.96 TeV, as measured by the DO 0 and CDL 0 
experiments at the Tevatron. Neither the ATLAS nor the CMS results show any statistically 
significant rapidity dependence of the cross section in the central region. The CMS analyses 
cover the px range up to 50 GeV, while the ATLAS results go to 70 GeV. 

In this Letter we present a measurement of the differential production cross sections of the three 
lowest-mass Y(nS) states in pp collisions at = 7 TeV up to px = 100 GeV, reaching higher 
Px than previous measurements. We measure the px dependence of the Y(nS) differential cross 
section times the branching fraction to using the 2011 data set, corresponding to an inte¬ 

grated luminosity of 4.9 fb~^. The measured cross sections include feeddown from higher bb 
excitations. 

Measurements of S-wave bb mesons provide an important probe of quantum chromod}mamics 
(QGD). There are several models that predict differential cross section shapes at high Y(nS) px 
in pp collisions. A common feature of all the models is that different contributing terms have 
different px variations, some of which are power-law forms. The nonrelativistic QGD (NRQGD) 
approach ||8l |9l uses an effective field theory to factorize the perturbative term and nonpertu- 
bative long-distance matrix element (LDME) terms. A good description of early LHG results 
for Y(1S) production for px < 30 GeV was achieved using NRQGD with next-to-leading-order 
(NLO) corrections [iHl. However, there are theoretical corrections to perturbative NRQGD that 
have characteristic power-law behavior at high px, and measurements at high px can help to 
clarify the theoretical picture ilTTlfT^ . The NLO NRQGD calculation has recently been extended 
to treat all three Y(nS) states 1131 . The updated calculation includes not only NLO terms but 
also uses LDMEs computed using only high-px data. Golor singlet models (GSM) with higher- 
order px-dependent corrections Iil4l and the /cx-fectorization model IH^ are consistent with 
data from the LHG for px approaching 50 GeV. A recent analysis of quarkonium polarization 
and production measurements found that raising px thresholds stabilizes the fits in evaluating 
the LDMEs |[T6] . At higher px different corrections become dominant in these models. New 
data at high px will challenge all the current approaches. 


2 CMS detector 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter having a 3.8 T field. Within the superconducting solenoid volume are a silicon pixel and 
strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator 
hadron calorimeter. Muons are measured in gas-ionization detectors embedded in the steel 
flux-return yoke outside the solenoid, with detection planes are made using three technolo¬ 
gies: drift tubes, cathode strip chambers, and resistive-plate chambers. Muons are measured in 
the pseudorapidity range |p| < 2.4. 

The silicon tracker measures charged particles within the pseudorapidity range \rj\ < 2.5. It 
consists of 1440 silicon pixel and 15148 silicon strip detector modules and provides a typical 
transverse impact parameter resolution of 25-90 pm. Matching muons to tracks measured in 
the silicon tracker results in a transverse momentum resolution between 1% and 2.8%, for px 
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values up to 100 GeV flT l. 

The first level of the CMS trigger system, composed of custom hardware processors, uses in¬ 
formation from the calorimeters and muon detectors to select the most interesting events in a 
fixed time interval of less than 4 ^s. The high-level trigger processor farm further decreases 
the event rate from around 100 kHz to around 400 Hz, before data storage. A more detailed 
description of the CMS detector, together with a definition of the coordinate system used and 
the relevant kinematic variables, can be found in Ref. flS l. 


3 Differential cross section measurement methodology 

Event selection starts with a dimuon trigger involving the silicon tracker and muon systems. 
The trigger, which is exposed to the full integrated luminosity, requires at least two muons 
with dimuon rapidity \y\ < 1.25, dimuon invariant mass 8.5 < < 11.5 GeV, and a dimuon 

vertex fit with a probability >0.5%. The trigger selects only pairs of muons that bend away 
from each other in the magnetic field ("seagull selection"), i.e., events for which the difference 
in azimuthal angle between the positively charged and negatively charged muons is less than 
zero. Requiring that muon trajectories do not cross in the transverse plane improves the muon 
efficiency. Trigger pj thresholds varied from 5-9 GeV as the beam conditions changed. Offline 
selection criteria required px > 10 GeV, \y\ < 1.2, and a dimuon vertex fit probability >1%. 
Standard GMS quality requirements are used to identify muons and muons are restricted to 
\y{y)\ < 1.6. The muon tracks are required to have at least ten hits in the silicon tracker, at 
least one hit in the silicon pixel detector, and be matched with at least one segment of the muon 
system. The muon track fit quality must have a per degree of freedom of less than 1.8. The 
distance of the track from the closest primary vertex must be less than 15 cm in the longitudinal 
direction and 3 cm in the transverse direction. The following kinematic requirements are also 
imposed to ensure accurate muon detection efficiency evaluation: 

PT(p)>3GeV for 1.4 < |p(p)| < 1.6, 

Pt(p) > 3.5GeV for 1.2 < |p(p)| < 1.4, (1) 

Pt(p) > 4.5GeV for |p(p)| < 1.2. 


The differential cross sections are measured for two rapidity ranges: \y\ < 0.6 and 0.6 < \y\ < 
1.2, as well as for the entire range \y\ < 1.2. In each rapidity range the data are binned in px, 
with bin edges at 2 GeV intervals between 10 and 40 GeV, then wider bins with edges at 43,46, 
50, 55, 60,70, and 100 GeV. 

The Y(nS) differential cross section times dimuon branching fraction, integrated over either of 
the two \y\ ranges and in a given px bin of width Apx, is 


dcr^pp ^ Y(nS)^ 


dpx 


B[Y{nS) —> p^p ) = 


^^Y{nS) 


(Pt) 


LApxe^;, (px)»4,(px)esgevp' 


( 2 ) 


I range 


where is the fitted number of Y (nS) events from the dimuon invariant mass distribution 

in a Px bin for the selected \y\ range, e^^ is the dimuon efficiency, L is the integrated luminosity, 
A is the polarization-corrected acceptance, Cgg is the efficiency of the seagull selection, and Cvp 
is the efficiency of the dimuon vertex x^ probability requirement. The efficiency and acceptance 
determinations are described below. 
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The total yield for the three Y(nS) states in the rapidity range \y\ <1.2 are 412900 ± 600 

Y(1S) events, 151700 ± 400 Y(2S) events, and 111 100 ± 300 Y(3S) events, where the uncer¬ 
tainties are statistical only The fine granularity of the CMS tracker kept the efficiency indepen¬ 
dent of changes in the LHC instantaneous luminosity throughout the ^/s = 7 TeV operations. 

3.1 Efficiency factors 

The dimuon efficiency for a given event is parameterized as: 

= ei[pT(Hi),?/(f/i)] e2[vj{n)>V{n)] (3) 

where is the overall single-muon quality and trigger efficiency The kinematic 

dependence of the p factor was determined in a study based on Monte Carlo (MC) simulation 
using EVTGEN |[T9]| with a detector simulation performed with Geant4 ||20|| . The parameter p 
accounts for the possibility that two genuine muons can be merged during the reconstruction 
or trigger selection, causing an inefficiency. It was found to depend on the quadrature sum 
of the differences Apj /(637 GeV), At], and 1.2A^ between the two muons. The MG simulation 
result was validated by measuring the p factor with Y(nS) events reconstructed using a data 
set that required only a single-muon trigger. In events such as those with pj < 50 GeV, where 
the muons are well separated, p = 1. For high-pj events of px > 80 GeV, where the muons are 
closer together, p drops to approximately 0.7. 

The single-muon efficiencies are measured using the tag-and-probe approach based on control 
samples in data, as described in Ref. IIITl . times the tracking efficiency (0.99 ± 0.01), determined 
from MG simulation. We assume that the dimuon efficiency within each Y (nS) mass region is 
the same for signal and background. The dimuon efficiency for a given (pj, |i/|) is obtained 
by averaging the calculated event dimuon efficiency for each data event in the bin. This 
is done separately for the three Y states, using a mass range of ±200 MeV for the Y(1S) and 
±100MeV for the higher-mass states. The narrower range for the Y(2S) and Y(3S) states is 
chosen because of the closeness in mass of these two states. The average efficiency, is 
typically 0.75-0.80. For all (px,!/) bins the systematic difference between averaging in or 
1 / £fiii can be neglected in comparison to the quoted systematic uncertainty due to the single 
muon efficiencies. To determine Ssg, we note that there is a 50% probability that a Y(nS) state 
will decay in the seagull configuration. It was verified in MG simulation that Cgg = 0.5. The 
efficiency Cvp for the dimuon vertex fit probability requirement is determined to be 0.99 ± 
0.01 from MG simulation, where the uncertainty is statistical. This efficiency was validated 
in data using events from a trigger that did not require vertex selection. We also computed 
the total acceptance and efficiency product in the MG simulation and compared it with the 
result based on the factorized approach. The results agreed over the entire px range of the 
measurement. 

3.2 Acceptance 

For each Y(nS) state the acceptance A is computed in each (px, \y\) bin and defined as the 
fraction of its dimuon decays that satisfy the single-muon kinematic selections given by Eq. Q. 
The acceptances are computed using generator-level muons, then repeated using reconstructed 
muons in the full simulation study. The results agree to better than 2% at all px values. Dif¬ 
ferences are contained within the systematic uncertainty band (Section 4.3) assigned for the 
muon reconstruction. To account for the effect of the Y(nS) polarization on the muon angular 
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distribution, each MC simulation event is weighted by an angular factor zv: 


w = 


471 \ 3 + Ae 


(l + Ag cos^ 0 + A^ sin^ 9 cos 2(j) + A^^ sin 26 cos cp), 


(4) 


where Ag, A^, Ag^ are the measured polarization parameters ||5l, 9 is the polar angle, and (p the 
azimuthal angle of the positively charged muon in the Y(nS) helicity frame (HX). The polar¬ 
ization was measured in the range 10 < px < 50 GeV in the same two rapidity bins as this 
analysis. The measured polarization parameters do not show a statistically significant depen¬ 
dence on Px- We linearly interpolate each of the measured polarization parameters in pj. Linear 
interpolation is also used for the 68.3% confidence level (CL) uncertainties in the polarization 
measurements to determine the uncertainty in the three parameters from the analysis. The po¬ 
larization parameters for px > 50 GeV are taken to be the average of the measured values for 
10 < Px < 50 GeV. The largest measured absolute uncertainty for each parameter is used for 
the extrapolated uncertainties because the spread in nominal values is small. The acceptance 
is computed initially using a flat px distribution within a bin, then reweighted after fitting the 
measured px distribution to a functional form (see Section 5). The acceptances in each px bin 
for the three rapidity intervals are given in Appendix]^ (Tables |A.7| - [A.15[ ) for the measured po¬ 
larization central value and the 68.3% CL uncertainties on the parameters |5]. In addition, we 
report the acceptance computed for the hypotheses of zero, 100% transverse, and 100% longitu¬ 
dinal polarization that correspond to the parameter values A^ = Ag^ = 0 and A^ = 0, +\, and — 1 
respectively. Because of the agreement in the acceptance when computed with generator-level 
and reconstructed muons, the cross section results reported here can be scaled to accommodate 
any other polarization by using a generator-level MC simulation with a given polarization. 


4 Yield determination procedure 


4.1 Lineshape determination 


The Y(nS) lineshape is determined using the measured muon momenta and their uncertainties, 
along with a generator-level simulated invariant mass (SIM) distribution including final-state 
radiation (FSR) effects. For events in a given (px, |y |) bin, the distribution of the dimuon invari¬ 
ant mass uncertainty ^ is computed from the muon track error matrices. 

In order to describe the Y(nS) SIM distribution without detector resolution effects, we simulate 
dimuon events for a given Y(nS) state using EVTGEN and compute the FSR using PHOTOS (221 
[23l . The standard PHOTOS minimum photon energy for the Y(nS) states is Ri50 MeV, which 
is of the same order as our dimuon invariant mass uncertainty. To improve the description, 
we extend the photon energy spectrum down to 2 MeV using a fit of the SIM distribution to 
the QED mner-bremsstrahlung formula 1123] . The systematic uncertainties of the soft photon 
approximation in PHOTOS compared to exact QED calculations are discussed in Ref. II23I . For 
the range of photon energies expected in Y (nS) decays the systematic uncertainty is negligible. 

In each rapidity range, the Y(nS) lineshape for a given px bin is expressed by a probability 
density function (PDF) for the signal dimuon mass M^^,. This function Cw, ^m) is the 

average of N values of the dimuon mass smeared with a resolution 


2 N j 

-^(M^^; Ca,, 5m) = ~1^ V® 

iV y/2nCwU 


/IclQ 


(5) 


Each Y(nS) state is handled in the same fashion. Values of m, and are selected by randomly 
sampling the radiative mass function and the ^ distribution for that (px, |i/|) bin. Two correction 
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factors are common to all three Y(nS) peaks in a given {pj, \y\) bin: a width scale factor c^, to 
correct for any ^ scale difference between data and the MC simulation, and a mass-shift ^m, 
to correct for any difference in px scale between data and the MC simulation. We sample 
N = 25 000 {rrii, ^i) points per px bin, stored in a histogram with 0.25 MeV bins to smooth 
the fluctuations and retain shape features. This histogram gives the normalized, resolution- 
smeared mass PDF for a given Y(nS) state in a particular (px, |i/|) bin. The procedure was 
validated in MC simulation by generating the lineshape using a subset of generated Y(1S) 
events, then fitting the rest of the events with that lineshape. The fitted number of events was 
consistent with the generated number. 

4.2 Fitting for yieids 

To determine the yields of the three states in each px and \y\ range requires a fit to the dimuon 
mass distribution in every (px, \y\) bin. The total PDF for describes the signal and back¬ 
ground contributions to the dimuon invariant mass distribution using a signal PDF as defined 
in Eq. ^ for each of the Y(nS) states, plus a background function. Four background functions 
are studied: an exponential and a Chebyshev series with maximum order of 0,1, or 2. 

We measure the yield by performing an extended maximum-likelihood fit using ROOPlT 12^ 
to determine the number of signal events associated with each normalized signal PDF. To allow 
cancellation of some common uncertainties in the muon acceptance and efficiency calculation 
in the measurement of the ratios of Y (2S) and Y(3S) differential cross sections to that of the 
Y(IS), we perform an additional fit normalized to the Y(IS) yield. For each px bin the optimal 
background function is determined using the Akaike Information Criterion (AIC) Il25ll . taking 
the function with the largest relative probability, as discussed in Ref. I 26l . This method is 
similar to a maximum-likelihood evaluation, but it adds a term equal to twice the number of 
free parameters in the fit, thus penalizing addition of free parameters. The parameters Cw and 
Sm are determined from the fit for each px bin. Typical values and corresponding uncertainties 
for Ciu and 3m are 1.04 ± 0.01 and 3 ± 1 MeV, respectively The fit correlation matrix shows that 
their influence on the yields is a small fraction of the statistical uncertainty in each yield. 

The plots in Fig. show two examples of fitting the dimuon invariant mass distribution using 
the lineshape method. The lower plots show the pull, (N^ata ~ /c^data/ in each dimuon mass 

bin, where N^ata is the observed number of events in the bin, is the integral of the fitted 
signal and background function in that bin, and the uncertainty u^ata is the Poisson statistical 
uncertainty. As can be seen in Fig. the lineshape description represents the data well, even at 
high Px and large rapidity. 

4.3 Systematic uncertainties 

The overall systematic uncertainty in the cross section for a given (px, \y\) bin includes uncer¬ 
tainties from the background fit method, the lineshape determination, the dimuon efficiency, 
the acceptance variations due to varying the polarization parameters within their 68.3% CF 
ranges, and the integrated luminosity. The systematic uncertainty from the background func¬ 
tion is estimated using the maximum difference in yields among background functions with 
an AIC probability above 5% Il25ll26l relative to the best background choice. An upper limit of 
1% on the systematic uncertainty from the lineshape function determination for all three Y(nS) 
states and all (px, \y\) bins is estimated by varying the width of the mass region in which the 
mass resolution parameter ^ is determined. The efficiency systematic is evaluated by modify¬ 
ing event by event, using the ±1 standard deviation values from the tag-and-probe mea¬ 
surements [SJ. There is a 1% systematic uncertainty to account for small variations in as 
a function of observed in the data. The measured p factor values from the experimental 
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4.9 fb'^ {7 TeV) 



4.9 fb'^ (7 TeV) 



Figure 1: Results of the fits to the dimuon invariant mass distribution for events in two bins: 
(a): \y\ < 0.6, 10 < px < 12GeV and (b): 0.6 < \y\ < 1.2, 50 < pj < 55GeV. The solid line 
is the result of the full fit. The dash-dotted line is the Y(1S) signal fit, the long-dashed line 
is the Y(2S) signal fit, and the dotted line is the Y(3S) signal fit. The short-dashed line is the 
background contribution. The lower plots show the pull for each mass bin. 


determination and from MG simulation agree over the full px range. We assign a systematic 
uncertainty for p of 0.5-5%, which equals the full difference between the MG simulation and 
the experimental measurement. We compute the acceptance systematic uncertainty by rais¬ 
ing and lowering all three polarization parameters by their interpolated 68.3% GL values from 
Ref. |5f. The resulting 5-8% change in the acceptance is used as the systematic uncertainty in 
the acceptance as tabulated in Appendix [^(Tables |A.7|jA.15| . The total systematic uncertainty 
is found from the quadrature sum of the individual systematic uncertainties. It is comparable 
to or smaller than the statistical uncertainty for px > 40 GeV. There is a 2.2% uncertainty 1271 
from the integrated luminosity determination that applies to all px bins. This uncertainty is not 
included in the uncertainties displayed in the figures or given in the tables. 


5 Results 

The measured Y (nS) differential cross sections versus px are shown in Fig.|^over the full rapid¬ 
ity range \y\ < 1.2. The vertical bars on the points in Fig. |^show the statistical and systematic 
uncertainties added in quadrature. Earlier GMS measurements 0 are shown for comparison, 
scaled by 0.5 to account for the smaller \y \ range in the latest measurement, where the scaling 
assumes that the rapidity distribution is flat. The Y (nS) differential cross sections peak near 
Px = 4 GeV, as seen in Fig. Their shape can be described by an exponential function for 
10 ^ Px ^ 20 GeV, while for px > 20 GeV the data lie above the exponential and the slope 
changes. Therefore, we fit the high-px measurements for each Y(nS) state using a power-law 
parametrization: 


dtr^pp —> Y(nS)^ 


dpx 


B(Y{nS) 


P P = 


I range 



( 6 ) 


where A is a normalization with units of pb/GeV. The value of po is fixed to 20 GeV and has 
no influence on the exponent a, which describes the curvature of the function. The differential 
cross section fits are evaluated using the integral value of the function over the px range of 
each bin, and the results are given in Table The bin centers are determined by the functional- 
weight method described in ||28l, using the exponential fit for px < 20 GeV and the power-law 
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form in Eq. for px > 20 GeV. Shifts from the pj-weighted mean values are negligible in all 
except the highest-px bin, where using the functional weight moves the bin center from 79 to 
82 GeV. Tables A.l - jA-Sl in Appendix [A| give the measured values shown in Fig. as well as for 
the two rapidity ranges \y\ < 0.6 and 0.6 < lyl < 1-2. 


To illustrate the quality of this functional description. Fig. (b) shows the fit results for the 
Y(1S) state with \y\ < 1.2. The solid line is the power-law fit for px > 20GeV. The dashed line 
is the exponential fit for 10 < px < 20 GeV. The lower plot shows, for each px bin, the pull 
determined from the differential cross section value in a (px, |i/|) bin and its total uncertainty. 


Table 1: The values of the parameters in Eq. (|^ from the power-law fit to Y(1S) events with 
Px > 20 GeV and \y\ <1.2, along with the value and the number of degrees of freedom 



Y(1S) 

Y(2S) 

Y(3S) 

A 

14.00 ± 0.75 

6.88 ± 0.48 

4.01 ± 0.30 

OL 

5.75 ± 0.07 

5.62 ± 0.10 

5.26 ± 0.10 

G 

0.45 ± 0.13 

0.62 ± 0.18 

0.26 ± 0.15 


8.7 

11 

15 

nd 

14 

14 

14 




Figure 2: (a) The Y (nS) differential px cross sections times dimuon branching fractions for \y\ < 
1.2. The Y (2S) and Y (3S) measurements are scaled by 0.1 and 0.01, respectively, for display pur¬ 
poses. The vertical bars show the total uncertainty, excluding the systematic uncertainty in the 
integrated luminosity. The horizontal bars show the bin widths. Previous GMS measurements 
for \y\ < 2.4 are shown as cross-hatched areas [IJ]- These results have been scaled by 0.5 to 
account for the smaller \y\ range in the latest measurement, where the scaling assumes that the 
rapidity distribution is flat. The solid lines are the NFO calculations from Ref. liT3l extended 
by the authors to cover the range px < 100 GeV. (b) Details of the parametrized cross section 
fit described in the text for Y(1S) with \y\ < 1.2. In this plot the solid line is the result of the 
power-law fit (see Eq. for px > 20 GeV. The dashed line shows an exponential fit to the data 
for 10 < Px < 20 GeV. The lower plot shows the pulls of the fit as defined in the text. 

Next, we consider the px dependence of the ratios of the Y (nS) production cross sections times 
their dimuon branching fractions. The yield fits are redone to compute explicitly the yield ratio 
r 2 i for Y(2S) to Y(IS) and for Y(3S) to Y(IS). The efficiency ratio is computed for each (px, 
\y\) bin. The polarization-weighted acceptance and its uncertainty is computed for each state 
separately, and the uncertainties are added in quadrature to determine the uncertainty in the 
ratio. The corrected yield ratios are R„i(px, \y\) = rniipj, \y\) (AiCi)/ (Ane„), where n = 2, 3. 
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6 Discussion 


The measured corrected ratios are shown in Fig.|^and given in Appendix [^(Tables |A.4| - [A.6| |. 
The rapid rise of both ratios for px < 20 GeV slows significantly for px > 20 GeV. The curves 
on the ratio plots are the ratios of the corresponding fitted functions from the individual Y(nS) 
differential cross section fits (exponential for px < 20 GeV, power-law for px > 20 GeV). The 
curves confirm that the change in ratios occurs in the same px range in which dc"/dpx also 
changes behavior. 

The measurements for the ratio R 31 in Fig. |^b), found in the supplementary material, can be 
fit to a linear function and to a constant in order to quantify the visual evidence that the Y(3S) 
production is harder than that of the Y(IS). The linear fit to measurements with px > 20 GeV 
has probability 0.22, while the fit to a constant has ^ probability 2.6 xlO^^. Thus, with 
relative probability 85000:1, we can say that Y(3S) production is harder than that of the Y(IS). 
The Y(2S) /Y(IS) production ratio versus px has a similar trend, but the statistical uncertainties 
are too large to make a definite statement. 




Figure 3: Measured differential cross section ratios as a function of px. Gorrected yield ratios: 
(a) R 21 ; (b) R 31 . The dashed line is the ratio of the exponential fits to the individual differential 
cross sections for 10 < px < 20 GeV. The solid line is the ratio of the corresponding power-law 
fits for Px > 20 GeV. 


6 Discussion 

Theoretical predictions for the Y(nS) differential cross sections have been previously compared 
to the first LHG cross section measurements lUOlH^IT^ . A more recent GMS measurement HJ] 
included the currently available predictions from the GSM JlTl, valid for px < 35 GeV, and an 
unpublished NRQGD prediction that covers the range px < 30 GeV. The NRQGD + NLO anal¬ 
ysis from Ref. lfT3]| describes Y(nS) production at Tevatron and LHG energies for px < 50 GeV. 
An extension of these predictions to px = 100 GeV is compared to the GMS measurements in 
Fig.|^(a). The calculations describe the trends of the data for all three Y(nS) states. 

The color evaporation model (GEM), a variant of the GSM, predicts that above a minimum 
Px ~ My(is)/ all bottomonium states should have the same px dependence 129 1. The measured 
ratios of the differential cross sections as a function of px in Fig. show that this is not the case 
for px less than about 40 GeV. 

Ghanging the Y (nS) px threshold for the data used in calculating the NRQGD predictions re¬ 
sults in different LDMEs ||T01I301Im1. Recent theoretical work [121 [161 has demonstrated the 
impact of varying the px thresholds in NRQGD analyses to study different production ampli¬ 
tude behavior. These new GMS data provide a significant extension of the px range that can 
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be used in evaluating matrix elements and studying pj-dependent corrections in NRQCD and 
other models. The new results on Y(3S) production are sufficiently accurate to allow one to 
focus model building of the px behavior on that state, for which feeddown contributions come 
only from the Xbi3P)- 


7 Summary 

Measurements of the differential production cross sections as a function of px for the Y(1S), 
Y(2S), and Y(3S) states in pp collisions at ^ys = 7TeV have been presented, based on a data 
sample corresponding to an integrated luminosity of 4.9 fb^^ collected by the CMS experiment 
at the LHC. Not only do these measurements significantly improve the precision of the results 
in previously analyzed px ranges ID-lSl, they also extend the maximum px range from 70 to 
100 GeV. Evidence has been presented for the first time of the power-law nature of the px 
distributions for all three Y(nS) states at high px. Combined with the CMS Y (nS) polarization 
results jSl, the new bottomonium measurements are a formidable challenge to our theoretical 
understanding of the production of heavy-quark bound states. 
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Table A.l: The px bin width, the weighted mean px within a bin, and the differential cross sec¬ 
tion times the dimuon branching fraction: {da/dpj)B{Y(nS) —> p^p^) for the Y(1S), Y(2S), 
and Y(3S) with 0 < \y\ < 0.6. The statistical uncertainty in the differential cross section is 
given as the percentage of the cross section in the format: astat/{da/ dpx) (%) and similarly 
for the systematic uncertainty. The percentage systematic uncertainty for a negative systematic 
shift is given in parentheses. The statistical uncertainties are derived from the fit to the dimuon 
mass spectrum. The systematic uncertainties are discussed in the text. The 2.2% systematic 
uncertainty in the integrated luminosity is not included. 





Y(1S) 



Y(2S) 



Y(3S) 


Px 

(Pt) 

— H 


f^syst 

— H 

^4tat 

^syst 

— B 

dpT^ 

f^stat 

f^syst 

da/dpj 

da/dpj 

da/dpj 

dtj /dpj 

da/dpj 

da/dpj 

GeV 

GeV 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

10-12 

11.0 

60936 

0.5 

6.8 (7.2) 

20036 

1.0 

8.6 (9.4) 

10951 

1.4 

10.0 (11.3) 

12-14 

13.0 

33828 

0.5 

5.1 (5.6) 

12072 

1.1 

6.0 (6.7) 

7181 

1.4 

7.1 (8.0) 

14-16 

15.0 

19670 

0.6 

4.8 (5.2) 

7621 

1.2 

5.2 (5.7) 

4596 

1.6 

5.7 (6.3) 

16-18 

17.0 

11504 

0.8 

4.6 (4.9) 

4858 

1.3 

5.1 (5.5) 

3097 

1.7 

5.6 (6.1) 

18-20 

19.0 

6914 

0.9 

4.5 (4.7) 

3034 

1.5 

4.9 (5.2) 

1966 

2.0 

5.2 (5.6) 

20-22 

21.0 

4235 

1.1 

4.2 (4.4) 

1932 

1.8 

4.6 (4.9) 

1423 

2.2 

5.0 (5.2) 

22-24 

23.0 

2721 

1.4 

4.1 (4.3) 

1203 

2.3 

4.5 (4.6) 

893 

2.7 

4.7 (4.9) 

24-26 

25.0 

1848 

1.7 

4.1 (4.2) 

849 

2.6 

4.6 (4.7) 

594 

3.3 

4.9 (5.1) 

26-28 

27.0 

1117 

2.2 

4.3 (4.4) 

537 

3.6 

5.7 (5.8) 

381 

4.4 

6.1 (6.2) 

28-30 

29.0 

845 

2.5 

4.4 (4.5) 

413 

4.0 

5.8 (5.8) 

314 

4.7 

6.0 (6.1) 

30-32 

31.0 

593 

2.8 

4.3 (4.4) 

288 

4.4 

5.0 (5.1) 

209 

5.2 

5.0 (5.2) 

32-34 

33.0 

420 

3.5 

4.7 (4.8) 

194 

6.0 

6.7 (6.8) 

162 

6.5 

6.4 (6.5) 

34-36 

35.0 

314 

3.8 

4.6 (4.6) 

158 

5.8 

5.9 (6.0) 

112 

7.1 

5.7 (5.8) 

36-38 

37.0 

209 

4.8 

4.5 (4.6) 

123 

6.5 

4.6 (4.7) 

92 

7.5 

4.6 (4.7) 

38-40 

39.0 

157 

5.4 

4.2 (4.2) 

86 

7.8 

4.4 (4.6) 

61 

9.5 

4.6 (4.7) 

40-43 

41.4 

114 

5.0 

4.0 (4.1) 

61 

7.6 

4.3 (4.4) 

42 

9.3 

4.4 (4.6) 

43-46 

44.4 

76 

6.3 

5.8 (5.8) 

36 

10.0 

9.5 (9.5) 

39 

9.4 

7.8 (7.9) 

46-50 

47.9 

49 

6.6 

4.2 (4.2) 

27 

9.6 

4.8 (4.9) 

22 

10.9 

5.1 (5.3) 

50-55 

52.4 

24 

9.2 

5.9 (5.8) 

17 

10.9 

4.8 (4.9) 

12 

13.6 

4.4 (4.6) 

55-60 

57.4 

15 

11.0 

6.3 (6.2) 

9.8 

14.8 

5.0 (5.1) 

8.5 

15.5 

4.6 (4.8) 

60-70 

64.6 

9.2 

10.1 

6.6 (6.5) 

4.7 

15.7 

7.3 (7.3) 

4.9 

14.8 

6.0 (6.1) 

70-100 

82.0 

2.3 

12.6 

11.9(11.8) 

1.2 

17.6 

8.9 (8.9) 

1.0 

19.1 

7.3 (7.2) 
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Table A.2: The pj bin width, the weighted mean pj within a bin, and the differential cross 


section times the dimuon branching fraction: {da/pj)B(Y{nS) —>• p^p 
and Y(3S) with 0.6 < |i/| < 1.2. The notation is the same as for Table 


forthe Y(1S),Y(2S), 


A.1 





Y(1S) 



Y(2S) 



Y(3S) 


Pj 

(Pt) 

dpT^ 

^stat 

^syst 

dpT^ 


^^syst 

dpT^ 

^^stat 

^syst 

da/dpj 

dcr/djjj 

dcr/dpj 

da/dpj 

dcr/dpj 

dcr/dpj 

GeV 

GeV 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

10-12 

11.0 

55260 

0.5 

6.1 (6.3) 

18371 

1.2 

7.8 (8.2) 

9855 

1.6 

8.9 (10.0) 

12-14 

13.0 

31331 

0.6 

4.1 (4.2) 

10973 

1.1 

4.9 (5.2) 

6741 

1.5 

5.5 (6.0) 

14-16 

15.0 

18063 

0.7 

3.8 (4.0) 

6685 

1.5 

4.6 (4.8) 

4298 

2.0 

5.5 (5.7) 

16-18 

17.0 

10481 

0.9 

3.8 (3.9) 

4105 

1.7 

4.9 (5.1) 

2759 

2.3 

4.8 (5.0) 

18-20 

19.0 

6286 

1.1 

3.7 (3.8) 

2624 

2.1 

4.6 (4.8) 

1871 

2.6 

4.6 (4.7) 

20-22 

21.0 

3875 

1.4 

3.8 (3.9) 

1746 

2.4 

4.7 (4.8) 

1220 

3.0 

5.0 (5.2) 

22-24 

23.0 

2439 

1.6 

3.6 (3.7) 

1188 

2.4 

3.8 (3.9) 

819 

3.0 

4.0 (4.2) 

24-26 

25.0 

1633 

1.9 

3.8 (3.9) 

865 

2.8 

4.0 (4.1) 

530 

3.7 

4.0 (4.2) 

26-28 

27.0 

1080 

2.2 

3.8 (3.9) 

551 

3.4 

4.5 (4.5) 

377 

4.3 

4.7 (4.8) 

28-30 

29.0 

765 

2.5 

3.7 (3.8) 

397 

3.9 

3.9 (4.0) 

274 

5.0 

4.2 (4.3) 

30-32 

31.0 

484 

3.8 

5.0 (5.1) 

229 

6.5 

9.0 (9.0) 

170 

7.8 

10.3 (10.3) 

32-34 

33.0 

371 

3.7 

6.0 (6.1) 

175 

6.0 

9.8 (9.8) 

123 

7.8 

10.3 (10.4) 

34-36 

35.0 

263 

4.3 

3.7 (3.8) 

151 

6.1 

4.0 (4.1) 

101 

8.1 

4.3 (4.4) 

36-38 

37.0 

193 

4.9 

3.6 (3.7) 

113 

7.0 

3.9 (4.0) 

84 

8.6 

4.3 (4.4) 

38-40 

39.0 

152 

5.7 

3.7 (3.8) 

87 

8.1 

4.0 (4.1) 

47 

11.9 

4.2 (4.3) 

40-43 

41.4 

98 

6.0 

5.0 (5.0) 

55 

8.4 

4.3 (4.4) 

43 

9.8 

4.2 (4.3) 

43-46 

44.4 

73 

6.5 

7.1 (7.1) 

37 

10.3 

9.9 (9.9) 

34 

11.1 

13.4 (13.4) 

46-50 

47.9 

51 

6.7 

4.0 (4.0) 

23 

11.0 

3.9 (4.0) 

16 

14.2 

4.0 (4.2) 

50-55 

52.4 

27 

8.6 

5.1 (5.1) 

14 

13.2 

4.5 (4.5) 

8.7 

18.3 

4.0 (4.2) 

55-60 

57.4 

15 

11.2 

5.1 (5.0) 

9.2 

15.5 

4.6 (4.6) 

11 

13.7 

4.5 (4.7) 

60-70 

64.6 

8.1 

10.8 

5.6 (5.6) 

3.6 

19.7 

5.1 (5.1) 

3.3 

20.9 

5.2 (5.3) 

70-100 

82.0 

1.6 

15.0 

7.4 (7.2) 

1.0 

24.0 

7.0 (6.9) 

1.3 

19.8 

6.6 (6.7) 
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Table A.3: The pj bin width, the weighted mean pj within a bin, and the differential cross 
section times the dimuon branching fraction: {da/pj)B{Y (nS) —> ) for the Y(IS), Y(2S), 

and Y(3S) with 0 < |y| < 1.2. The notation is the same as for Table 


A.l 





Y(1S) 



Y(2S) 



Y(3S) 


Pj 

(Pt) 

dpT^ 

^stat 

^syst 

dpT^ 


^^syst 

dpT^ 

^^stat 

^^syst 

da/dpj 

dcr/djjj 

dcr/dpj 

da/dpj 

dcr/dpj 

dcr/dpj 

GeV 

GeV 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

(fb/GeV) 

(%) 

(%) 

10-12 

11.0 

116415 

0.4 

6.4 (6.7) 

38540 

0.7 

8.1 (8.7) 

20882 

1.1 

9.3 (10.5) 

12-14 

13.0 

65266 

0.4 

4.6 (4.9) 

23088 

0.8 

5.4 (5.9) 

13947 

1.0 

6.2 (6.9) 

14-16 

15.0 

37778 

0.5 

4.3 (4.6) 

14321 

0.9 

4.7 (5.1) 

8909 

1.3 

5.2 (5.7) 

16-18 

17.0 

22008 

0.6 

4.1 (4.4) 

8969 

1.1 

4.7 (5.0) 

5873 

1.4 

5.0 (5.3) 

18-20 

19.0 

13212 

0.7 

4.1 (4.2) 

5665 

1.3 

4.6 (4.8) 

3842 

1.6 

4.7 (5.0) 

20-22 

21.0 

8116 

0.9 

4.0 (4.1) 

3683 

1.5 

4.4 (4.6) 

2648 

1.8 

4.7 (4.9) 

22-24 

23.0 

5162 

1.0 

3.9 (4.0) 

2393 

1.7 

4.1 (4.2) 

1713 

2.0 

4.3 (4.5) 

24-26 

25.0 

3483 

1.3 

3.9 (4.0) 

1715 

1.9 

4.2 (4.3) 

1124 

2.4 

4.4 (4.5) 

26-28 

27.0 

2197 

1.5 

4.0 (4.0) 

1089 

2.5 

4.6 (4.7) 

759 

3.1 

4.9 (5.0) 

28-30 

29.0 

1611 

1.8 

4.0 (4.1) 

811 

2.8 

4.6 (4.7) 

588 

3.4 

4.8 (5.0) 

30-32 

31.0 

1077 

2.3 

4.2 (4.3) 

517 

3.8 

5.7 (5.7) 

380 

4.5 

6.2 (6.3) 

32-34 

33.0 

791 

2.5 

4.7 (4.7) 

369 

4.2 

6.5 (6.6) 

286 

5.0 

6.5 (6.6) 

34-36 

35.0 

577 

2.8 

4.1 (4.1) 

309 

4.2 

4.7 (4.7) 

213 

5.3 

4.8 (4.9) 

36-38 

37.0 

402 

3.4 

4.0 (4.1) 

236 

4.8 

4.2 (4.3) 

176 

5.7 

4.4 (4.6) 

38-40 

39.0 

308 

3.9 

3.9 (3.9) 

173 

5.6 

4.2 (4.3) 

109 

7.5 

4.4 (4.5) 

40-43 

41.4 

212 

3.9 

4.1 (4.2) 

116 

5.6 

4.2 (4.3) 

86 

6.8 

4.3 (4.5) 

43-46 

44.4 

148 

4.5 

5.3 (5.3) 

73 

7.2 

7.5 (7.5) 

73 

7.2 

8.1 (8.2) 

46-50 

47.9 

100 

4.7 

4.0 (4.0) 

50 

7.3 

4.2 (4.3) 

38 

8.7 

4.5 (4.6) 

50-55 

52.4 

51 

6.3 

5.0 (5.0) 

31 

8.4 

4.5 (4.6) 

21 

11.0 

4.2 (4.4) 

55-60 

57.4 

29 

7.9 

5.4 (5.4) 

19 

10.7 

4.8 (4.8) 

20 

10.3 

4.6 (4.7) 

60-70 

64.6 

17 

7.4 

6.0 (5.9) 

8.3 

12.4 

6.1 (6.1) 

8.1 

12.2 

5.5 (5.6) 

70-100 

82.0 

3.9 

9.7 

8.9 (8.8) 

2.2 

14.6 

7.6 (7.5) 

2.3 

14.0 

6.9 (6.9) 
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Table A.4: The px bin width and corrected yield ratios R 21 arid R 31 for \y\ < 0.6. The sta¬ 
tistical uncertainty in the corrected yield ratios is given as the percentage of the ratio in the 
format: Cstat/Rni (%) {n = 2,3) and similarly for the systematic uncertainty. The statistical un¬ 
certainties are derived from the fit to the dimuon mass spectrum. The systematic uncertainties 
are discussed in the text. The 2.2% systematic uncertainty in the integrated luminosity is not 
included. 


Pt 

R21 

R21 

tTsyst 

R21 

R31 

R 31 

O'syst 

R31 

GeV 


(%) 

(%) 


(%) 

(%) 

10- 

-12 

0.33 

1.0 

7.2 

0.18 

1.4 

8.6 

12- 

-14 

0.36 

1.1 

5.9 

0.21 

1.5 

6.7 

14- 

-16 

0.39 

1.3 

4.9 

0.23 

1.7 

5.5 

16- 

-18 

0.42 

1.4 

4.7 

0.27 

1.8 

5.2 

18- 

-20 

0.44 

1.7 

4.1 

0.28 

2.1 

4.6 

20- 

-22 

0.46 

2.1 

3.8 

0.34 

2.4 

4.2 

22- 

-24 

0.44 

2.6 

3.4 

0.33 

3.0 

3.8 

24- 

-26 

0.46 

3.0 

3.4 

0.32 

3.6 

3.8 

26- 

-28 

0.48 

3.9 

3.9 

0.34 

4.6 

4.5 

28- 

-30 

0.49 

4.4 

3.8 

0.37 

5.0 

4.1 

30- 

-32 

0.49 

5.1 

3.5 

0.35 

5.8 

3.6 

32- 

-34 

0.46 

6.4 

4.1 

0.39 

6.9 

4.1 

34- 

-36 

0.50 

6.8 

3.6 

0.36 

7.8 

3.7 

36- 

-38 

0.59 

7.8 

3.3 

0.44 

8.6 

4.0 

38-40 

0.55 

9.2 

3.1 

0.39 

10.6 

3.5 

40-43 

0.54 

8.8 

3.0 

0.37 

10.3 

3.3 

43-46 

0.48 

11.4 

5.3 

0.51 

10.9 

3.9 

46- 

-50 

0.55 

11.3 

3.1 

0.45 

12.3 

3.5 

50- 

-55 

0.72 

13.4 

3.2 

0.51 

15.6 

4.6 

55- 

-60 

0.66 

17.4 

2.9 

0.57 

18.2 

4.8 

60- 

-70 

0.51 

17.8 

2.5 

0.53 

16.9 

2.8 

70- 

100 

0.53 

20.3 

4.1 

0.44 

21.9 

9.1 
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Table A.5: The pj bin width and corrected yield ratios R 21 arid R 31 for 0.6 < \y\ < 1.2. The 
notation is the same as for Table 


A.4 


Pt 

R 21 

R21 

f^syst 

R21 

R 31 

i^stat 

R31 

f^syst 

^31 

GeV 


(%) 

(%) 


(%) 

(%) 

10- 

-12 

0.33 

1.1 

5.5 

0.18 

1.6 

6.4 

12- 

-14 

0.35 

1.2 

3.9 

0.22 

1.5 

4.3 

14- 

-16 

0.37 

1.4 

3.6 

0.24 

1.9 

4.1 

16- 

-18 

0.39 

1.7 

3.4 

0.26 

2.2 

4.1 

18- 

-20 

0.42 

2.1 

3.4 

0.30 

2.5 

3.6 

20- 

-22 

0.45 

2.4 

3.2 

0.31 

3.0 

3.8 

22- 

-24 

0.49 

2.8 

2.8 

0.34 

3.3 

3.0 

24- 

-26 

0.53 

3.2 

2.8 

0.32 

4.0 

3.3 

26- 

-28 

0.51 

3.9 

3.1 

0.35 

4.7 

3.2 

28- 

0 

CO 

1 

0.52 

4.6 

2.9 

0.36 

5.5 

3.4 

0 

CO 

-32 

0.47 

6.4 

5.3 

0.35 

7.6 

6.6 

32- 

-34 

0.47 

6.8 

5.3 

0.33 

8.2 

5.9 

34- 

-36 

0.58 

7.3 

3.0 

0.39 

8.9 

3.5 

36- 

00 

CO 

1 

0.58 

8.2 

2.9 

0.44 

9.5 

3.4 

38-40 

0.57 

9.4 

2.9 

0.31 

12.8 

3.2 

40-43 

0.56 

9.8 

3.0 

0.44 

11.0 

4.5 

43-46 

0.51 

12.9 

6.8 

0.47 

15.5 

5.4 

46- 

-50 

0.45 

12.3 

2.6 

0.31 

15.0 

3.3 

50- 

-55 

0.51 

14.9 

2.6 

0.33 

19.4 

4.2 

55- 

-60 

0.63 

18.4 

4.6 

0.77 

16.6 

2.9 

60- 

-70 

0.46 

21.0 

4.7 

0.38 

23.2 

6.4 

70- 

100 

0.66 

25.8 

4.2 

0.75 

23.4 

9.7 
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Table A. 6 : The px bin width and corrected yield ratios R 21 and R 31 for \y\ < 1.2. The notation 
is the same as for Table 


A.4 


Pt 

GeV 

R21 

R21 

(%) 

f^syst 

R21 

(%) 

R31 

R31 

(%) 

f^syst 

^31 

(%) 

10-12 

0.33 

0.7 

6.2 

0.18 

1.0 

7.4 

12-14 

0.35 

0.8 

4.8 

0.21 

1.1 

5.4 

14-16 

0.38 

1.0 

4.1 

0.24 

1.3 

4.5 

16-18 

0.41 

1.1 

4.0 

0.27 

1.4 

4.3 

18-20 

0.43 

1.3 

3.6 

0.29 

1.6 

3.9 

20-22 

0.45 

1.6 

3.4 

0.33 

1.9 

3.7 

22-24 

0.46 

1.9 

3.0 

0.33 

2.2 

3.3 

24-26 

0.49 

2.2 

3.0 

0.32 

2.7 

3.4 

26-28 

0.50 

2.8 

3.2 

0.34 

3.3 

3.4 

28-30 

0.50 

3.2 

3.1 

0.36 

3.7 

3.5 

30-32 

0.48 

4.0 

3.3 

0.35 

4.6 

3.5 

32-34 

0.47 

4.7 

3.8 

0.36 

5.3 

4.0 

34-36 

0.53 

5.0 

3.1 

0.37 

5.9 

3.5 

36-38 

0.58 

5.7 

3.1 

0.44 

6.4 

3.4 

38-40 

0.56 

6.6 

3.0 

0.35 

8.2 

3.3 

40-43 

0.55 

6.6 

2.9 

0.40 

7.5 

3.2 

43-46 

0.49 

8.6 

4.5 

0.50 

8.9 

3.7 

46-50 

0.49 

8.4 

2.7 

0.37 

9.7 

3.2 

50-55 

0.59 

10.1 

2.6 

0.40 

12.5 

3.6 

55-60 

0.65 

12.6 

2.8 

0.65 

12.4 

3.0 

60-70 

0.49 

13.6 

2.5 

0.45 

13.9 

2.6 

70-100 

0.56 

16.0 

3.3 

0.50 

16.5 

7.0 
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Table A. 7: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(IS) state and \y\ < 0.6. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.31 

0.29 

0.32 

0.32 

0.26 

0.45 

12-14 

0.38 

0.36 

0.40 

0.39 

0.32 

0.54 

14-16 

0.44 

0.42 

0.45 

0.45 

0.36 

0.61 

16-18 

0.49 

0.48 

0.51 

0.50 

0.41 

0.67 

18-20 

0.53 

0.52 

0.55 

0.54 

0.45 

0.71 

20-22 

0.57 

0.56 

0.58 

0.57 

0.48 

0.74 

22-24 

0.60 

0.59 

0.61 

0.60 

0.51 

0.77 

24-26 

0.62 

0.61 

0.64 

0.62 

0.54 

0.79 

26-28 

0.64 

0.63 

0.66 

0.65 

0.56 

0.81 

28-30 

0.66 

0.64 

0.67 

0.66 

0.58 

0.83 

30-32 

0.67 

0.66 

0.69 

0.68 

0.60 

0.84 

32-34 

0.69 

0.67 

0.70 

0.69 

0.62 

0.85 

34-36 

0.70 

0.69 

0.71 

0.71 

0.63 

0.86 

36-38 

0.71 

0.70 

0.72 

0.72 

0.64 

0.87 

38-40 

0.72 

0.71 

0.74 

0.73 

0.66 

0.88 

40-43 

0.73 

0.72 

0.75 

0.74 

0.67 

0.89 

43-46 

0.75 

0.73 

0.76 

0.76 

0.69 

0.89 

46-50 

0.76 

0.75 

0.78 

0.77 

0.71 

0.90 

50-55 

0.78 

0.77 

0.80 

0.79 

0.73 

0.91 

55-60 

0.80 

0.78 

0.81 

0.80 

0.75 

0.92 

60-70 

0.82 

0.81 

0.83 

0.82 

0.77 

0.93 

70-100 

0.87 

0.86 

0.88 

0.87 

0.83 

0.96 
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Table A.8: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(IS) state and 0.6 < \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.29 

0.28 

0.30 

0.29 

0.24 

0.38 

12-14 

0.35 

0.34 

0.36 

0.35 

0.30 

0.46 

14-16 

0.42 

0.41 

0.43 

0.42 

0.36 

0.54 

16-18 

0.47 

0.46 

0.48 

0.47 

0.40 

0.60 

18-20 

0.52 

0.51 

0.53 

0.51 

0.44 

0.66 

20-22 

0.55 

0.54 

0.56 

0.55 

0.48 

0.70 

22-24 

0.59 

0.58 

0.60 

0.59 

0.51 

0.74 

24-26 

0.61 

0.60 

0.63 

0.62 

0.54 

0.77 

26-28 

0.64 

0.63 

0.65 

0.64 

0.57 

0.79 

28-30 

0.66 

0.64 

0.67 

0.66 

0.59 

0.81 

30-32 

0.67 

0.66 

0.69 

0.68 

0.61 

0.83 

32-34 

0.69 

0.67 

0.70 

0.70 

0.63 

0.84 

34-36 

0.70 

0.69 

0.71 

0.71 

0.64 

0.86 

36-38 

0.71 

0.70 

0.72 

0.73 

0.66 

0.87 

38-40 

0.72 

0.71 

0.74 

0.74 

0.67 

0.88 

40-43 

0.74 

0.73 

0.75 

0.75 

0.68 

0.88 

43-46 

0.75 

0.74 

0.77 

0.77 

0.70 

0.89 

46-50 

0.77 

0.76 

0.78 

0.78 

0.72 

0.90 

50-55 

0.79 

0.78 

0.80 

0.80 

0.74 

0.91 

55-60 

0.81 

0.80 

0.82 

0.81 

0.76 

0.92 

60-70 

0.83 

0.82 

0.84 

0.83 

0.78 

0.93 

70-100 

0.88 

0.87 

0.89 

0.88 

0.84 

0.96 
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Table A.9: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(IS) state and \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.30 

0.28 

0.31 

0.30 

0.25 

0.42 

12-14 

0.37 

0.35 

0.38 

0.37 

0.31 

0.50 

14-16 

0.43 

0.42 

0.44 

0.43 

0.36 

0.57 

16-18 

0.48 

0.47 

0.49 

0.48 

0.41 

0.64 

18-20 

0.52 

0.51 

0.54 

0.53 

0.45 

0.68 

20-22 

0.56 

0.55 

0.57 

0.56 

0.48 

0.72 

22-24 

0.59 

0.58 

0.61 

0.59 

0.51 

0.75 

24-26 

0.62 

0.61 

0.63 

0.62 

0.54 

0.78 

26-28 

0.64 

0.63 

0.65 

0.64 

0.56 

0.80 

28-30 

0.66 

0.64 

0.67 

0.66 

0.58 

0.82 

30-32 

0.67 

0.66 

0.69 

0.68 

0.60 

0.84 

32-34 

0.69 

0.67 

0.70 

0.70 

0.62 

0.85 

34-36 

0.70 

0.69 

0.71 

0.71 

0.64 

0.86 

36-38 

0.71 

0.70 

0.72 

0.72 

0.65 

0.87 

38-40 

0.72 

0.71 

0.74 

0.73 

0.66 

0.88 

40-43 

0.74 

0.72 

0.75 

0.75 

0.68 

0.88 

43-46 

0.75 

0.74 

0.77 

0.76 

0.70 

0.89 

46-50 

0.77 

0.76 

0.78 

0.78 

0.71 

0.90 

50-55 

0.79 

0.77 

0.80 

0.79 

0.73 

0.91 

55-60 

0.80 

0.79 

0.81 

0.81 

0.75 

0.92 

60-70 

0.82 

0.81 

0.84 

0.83 

0.78 

0.93 

70-100 

0.87 

0.86 

0.89 

0.88 

0.83 

0.96 
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Table A.10: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(2S) state and \y\ < 0.6. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.31 

0.30 

0.33 

0.34 

0.27 

0.48 

12-14 

0.37 

0.36 

0.39 

0.41 

0.33 

0.56 

14-16 

0.43 

0.41 

0.45 

0.47 

0.38 

0.63 

16-18 

0.48 

0.46 

0.50 

0.51 

0.42 

0.68 

18-20 

0.52 

0.51 

0.54 

0.55 

0.46 

0.73 

20-22 

0.56 

0.54 

0.57 

0.58 

0.49 

0.76 

22-24 

0.59 

0.57 

0.61 

0.61 

0.52 

0.79 

24-26 

0.61 

0.59 

0.63 

0.63 

0.55 

0.81 

26-28 

0.64 

0.62 

0.65 

0.66 

0.58 

0.83 

28-30 

0.65 

0.63 

0.67 

0.67 

0.59 

0.84 

30-32 

0.66 

0.64 

0.68 

0.69 

0.61 

0.85 

32-34 

0.68 

0.66 

0.70 

0.71 

0.63 

0.87 

34-36 

0.69 

0.67 

0.71 

0.72 

0.64 

0.88 

36-38 

0.70 

0.68 

0.72 

0.73 

0.65 

0.88 

38-40 

0.71 

0.69 

0.73 

0.74 

0.66 

0.89 

40-43 

0.73 

0.71 

0.75 

0.76 

0.69 

0.90 

43-46 

0.74 

0.72 

0.76 

0.77 

0.70 

0.91 

46-50 

0.76 

0.74 

0.78 

0.79 

0.72 

0.92 

50-55 

0.77 

0.76 

0.79 

0.80 

0.74 

0.92 

55-60 

0.79 

0.78 

0.81 

0.82 

0.76 

0.93 

60-70 

0.82 

0.80 

0.83 

0.84 

0.78 

0.95 

70-100 

0.87 

0.85 

0.88 

0.88 

0.84 

0.97 
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Table A.11: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(2S) state and 0.6 < \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A(cr ) 

A(unpol) 

A(T) 

A(L) 

10-12 

0.28 

0.27 

0.30 

0.30 

0.25 

0.39 

12-14 

0.35 

0.34 

0.37 

0.36 

0.31 

0.47 

14-16 

0.42 

0.41 

0.43 

0.42 

0.36 

0.54 

16-18 

0.47 

0.46 

0.49 

0.47 

0.41 

0.60 

18-20 

0.51 

0.50 

0.53 

0.52 

0.45 

0.66 

20-22 

0.54 

0.53 

0.56 

0.56 

0.49 

0.70 

22-24 

0.56 

0.55 

0.58 

0.59 

0.52 

0.74 

24-26 

0.58 

0.57 

0.60 

0.62 

0.54 

0.77 

26-28 

0.61 

0.60 

0.63 

0.65 

0.57 

0.79 

28-30 

0.64 

0.62 

0.66 

0.67 

0.60 

0.82 

30-32 

0.66 

0.64 

0.68 

0.69 

0.61 

0.84 

32-34 

0.68 

0.66 

0.69 

0.70 

0.63 

0.85 

34-36 

0.69 

0.68 

0.71 

0.72 

0.65 

0.87 

36-38 

0.72 

0.70 

0.73 

0.74 

0.67 

0.88 

38-40 

0.72 

0.71 

0.74 

0.75 

0.68 

0.89 

40-43 

0.74 

0.72 

0.76 

0.76 

0.69 

0.89 

43-46 

0.76 

0.74 

0.78 

0.78 

0.71 

0.91 

46-50 

0.77 

0.75 

0.79 

0.79 

0.73 

0.92 

50-55 

0.79 

0.78 

0.81 

0.81 

0.75 

0.93 

55-60 

0.81 

0.80 

0.83 

0.83 

0.77 

0.94 

60-70 

0.83 

0.82 

0.85 

0.84 

0.79 

0.95 

70-100 

0.88 

0.87 

0.90 

0.89 

0.85 

0.97 
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Table A.12: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(2S) state and \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.30 

0.28 

0.32 

0.32 

0.26 

0.44 

12-14 

0.36 

0.35 

0.38 

0.38 

0.32 

0.52 

14-16 

0.42 

0.41 

0.44 

0.44 

0.37 

0.59 

16-18 

0.47 

0.46 

0.49 

0.49 

0.41 

0.64 

18-20 

0.52 

0.50 

0.53 

0.53 

0.46 

0.69 

20-22 

0.55 

0.53 

0.56 

0.57 

0.49 

0.73 

22-24 

0.58 

0.56 

0.59 

0.60 

0.52 

0.76 

24-26 

0.60 

0.58 

0.61 

0.62 

0.54 

0.79 

26-28 

0.62 

0.61 

0.64 

0.65 

0.57 

0.81 

28-30 

0.64 

0.63 

0.66 

0.67 

0.59 

0.83 

30-32 

0.66 

0.64 

0.68 

0.69 

0.61 

0.85 

32-34 

0.68 

0.66 

0.70 

0.71 

0.63 

0.86 

34-36 

0.69 

0.67 

0.71 

0.72 

0.65 

0.87 

36-38 

0.71 

0.69 

0.72 

0.73 

0.66 

0.88 

38-40 

0.71 

0.70 

0.73 

0.74 

0.67 

0.89 

40-43 

0.73 

0.72 

0.75 

0.76 

0.69 

0.90 

43-46 

0.75 

0.73 

0.77 

0.77 

0.71 

0.91 

46-50 

0.76 

0.75 

0.78 

0.79 

0.72 

0.92 

50-55 

0.78 

0.77 

0.80 

0.80 

0.74 

0.93 

55-60 

0.80 

0.79 

0.82 

0.82 

0.76 

0.93 

60-70 

0.82 

0.81 

0.84 

0.84 

0.79 

0.95 

70-100 

0.87 

0.86 

0.89 

0.88 

0.84 

0.97 
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Table A.13: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(3S) state and \y\ < 0.6. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.34 

0.32 

0.37 

0.35 

0.28 

0.49 

12-14 

0.39 

0.37 

0.42 

0.42 

0.34 

0.58 

14-16 

0.44 

0.42 

0.46 

0.47 

0.39 

0.64 

16-18 

0.49 

0.47 

0.51 

0.52 

0.43 

0.69 

18-20 

0.52 

0.50 

0.55 

0.55 

0.46 

0.73 

20-22 

0.55 

0.54 

0.58 

0.59 

0.50 

0.77 

22-24 

0.59 

0.57 

0.60 

0.62 

0.53 

0.79 

24-26 

0.61 

0.59 

0.63 

0.64 

0.55 

0.82 

26-28 

0.63 

0.62 

0.66 

0.67 

0.58 

0.84 

28-30 

0.65 

0.63 

0.67 

0.68 

0.60 

0.85 

30-32 

0.67 

0.65 

0.69 

0.70 

0.61 

0.86 

32-34 

0.69 

0.67 

0.71 

0.71 

0.63 

0.87 

34-36 

0.70 

0.68 

0.72 

0.73 

0.65 

0.88 

36-38 

0.72 

0.70 

0.74 

0.74 

0.67 

0.89 

38-40 

0.73 

0.70 

0.75 

0.75 

0.67 

0.90 

40-43 

0.74 

0.72 

0.76 

0.76 

0.69 

0.91 

43-46 

0.75 

0.73 

0.78 

0.78 

0.71 

0.91 

46-50 

0.77 

0.75 

0.79 

0.79 

0.72 

0.92 

50-55 

0.79 

0.77 

0.81 

0.81 

0.74 

0.93 

55-60 

0.81 

0.79 

0.83 

0.82 

0.77 

0.94 

60-70 

0.82 

0.81 

0.85 

0.84 

0.79 

0.95 

70-100 

0.87 

0.85 

0.89 

0.88 

0.84 

0.98 
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Table A.14: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(3S) state and 0.6 < \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.31 

0.29 

0.33 

0.30 

0.25 

0.40 

12-14 

0.35 

0.34 

0.37 

0.36 

0.31 

0.47 

14-16 

0.40 

0.39 

0.41 

0.42 

0.36 

0.54 

16-18 

0.44 

0.43 

0.46 

0.48 

0.41 

0.60 

18-20 

0.49 

0.48 

0.50 

0.52 

0.46 

0.66 

20-22 

0.53 

0.52 

0.55 

0.56 

0.49 

0.70 

22-24 

0.58 

0.56 

0.59 

0.59 

0.52 

0.73 

24-26 

0.61 

0.60 

0.63 

0.62 

0.55 

0.77 

26-28 

0.63 

0.61 

0.65 

0.64 

0.57 

0.79 

28-30 

0.65 

0.63 

0.67 

0.67 

0.60 

0.82 

30-32 

0.66 

0.65 

0.68 

0.69 

0.62 

0.84 

32-34 

0.68 

0.66 

0.70 

0.71 

0.64 

0.85 

34-36 

0.69 

0.67 

0.71 

0.72 

0.65 

0.87 

36-38 

0.70 

0.68 

0.72 

0.74 

0.67 

0.88 

38-40 

0.72 

0.70 

0.74 

0.76 

0.69 

0.89 

40-43 

0.73 

0.72 

0.76 

0.77 

0.70 

0.90 

43-46 

0.75 

0.73 

0.77 

0.78 

0.72 

0.91 

46-50 

0.77 

0.75 

0.79 

0.80 

0.73 

0.92 

50-55 

0.79 

0.77 

0.81 

0.81 

0.75 

0.93 

55-60 

0.81 

0.79 

0.83 

0.83 

0.77 

0.94 

60-70 

0.83 

0.82 

0.85 

0.85 

0.80 

0.95 

70-100 

0.88 

0.87 

0.90 

0.89 

0.85 

0.98 
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Table A.15: The dimuon acceptance A calculated using the CMS measured polarization and 
its positive and negative uncertainties A{cr^), and the values of the acceptance assuming no 
polarization A(unpol), transverse polarization A{T), and longitudinal polarization A{L) for 
the Y(3S) state and \y\ < 1.2. 


PT [GeV] 

A 

A{a+) 

A{(t ) 

A(unpol) 

A{T) 

A{L) 

10-12 

0.32 

0.30 

0.35 

0.33 

0.27 

0.44 

12-14 

0.37 

0.36 

0.39 

0.39 

0.33 

0.52 

14-16 

0.42 

0.40 

0.44 

0.45 

0.37 

0.59 

16-18 

0.47 

0.45 

0.48 

0.50 

0.42 

0.65 

18-20 

0.51 

0.49 

0.52 

0.54 

0.46 

0.69 

20-22 

0.54 

0.53 

0.56 

0.57 

0.49 

0.73 

22-24 

0.58 

0.57 

0.60 

0.60 

0.53 

0.76 

24-26 

0.61 

0.59 

0.63 

0.63 

0.55 

0.79 

26-28 

0.63 

0.61 

0.65 

0.65 

0.58 

0.81 

28-30 

0.65 

0.63 

0.67 

0.68 

0.60 

0.83 

30-32 

0.67 

0.65 

0.69 

0.69 

0.62 

0.85 

32-34 

0.68 

0.66 

0.70 

0.71 

0.63 

0.86 

34-36 

0.69 

0.67 

0.72 

0.72 

0.65 

0.87 

36-38 

0.71 

0.69 

0.73 

0.74 

0.67 

0.89 

38-40 

0.72 

0.70 

0.75 

0.75 

0.68 

0.89 

40-43 

0.74 

0.72 

0.76 

0.77 

0.70 

0.90 

43-46 

0.75 

0.73 

0.77 

0.78 

0.71 

0.91 

46-50 

0.77 

0.75 

0.79 

0.79 

0.73 

0.92 

50-55 

0.79 

0.77 

0.81 

0.81 

0.75 

0.93 

55-60 

0.81 

0.79 

0.83 

0.83 

0.77 

0.94 

60-70 

0.83 

0.81 

0.85 

0.85 

0.79 

0.95 

70-100 

0.88 

0.86 

0.90 

0.89 

0.84 

0.98 
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